% Plot the data in
function[] = imagesc_on_map(CUr,CUi,CVr,CVi)

theta = [0.125 : 0.25 : 180.]';
phi   = [0.125 : 0.25 : 360.]';

map = [    0.2510         0    0.2510
    0.2331         0    0.3045
    0.2151         0    0.3580
    0.1972         0    0.4115
    0.1793         0    0.4650
    0.1613         0    0.5185
    0.1434         0    0.5720
    0.1255         0    0.6255
    0.1076         0    0.6790
    0.0896         0    0.7325
    0.0717         0    0.7860
    0.0538         0    0.8395
    0.0359         0    0.8930
    0.0179         0    0.9465
         0         0    1.0000
    0.0091    0.0909    1.0000
    0.0182    0.1818    1.0000
    0.0273    0.2727    1.0000
    0.0364    0.3636    1.0000
    0.0455    0.4545    1.0000
    0.0545    0.5455    1.0000
    0.0636    0.6364    1.0000
    0.0727    0.7273    1.0000
    0.0818    0.8182    1.0000
    0.0909    0.9091    1.0000
    0.1000    1.0000    1.0000
    0.2500    1.0000    1.0000
    0.4000    1.0000    1.0000
    0.5500    1.0000    1.0000
    0.7000    1.0000    1.0000
    0.8500    1.0000    1.0000
    1.0000    1.0000    1.0000
    1.0000    1.0000    0.8333
    1.0000    1.0000    0.6667
    1.0000    1.0000    0.5000
    1.0000    1.0000    0.3333
    1.0000    1.0000    0.1667
    1.0000    1.0000         0
    1.0000    0.9333         0
    1.0000    0.8667         0
    1.0000    0.8000         0
    1.0000    0.7333         0
    1.0000    0.6667         0
    1.0000    0.6000         0
    1.0000    0.5333         0
    1.0000    0.4667         0
    1.0000    0.4000         0
    1.0000    0.3333         0
    1.0000    0.2667         0
    1.0000    0.2000         0
    1.0000    0.1333         0
    1.0000    0.0667         0
    1.0000         0         0
    1.0000         0    0.0909
    1.0000         0    0.1818
    1.0000         0    0.2727
    1.0000         0    0.3636
    1.0000         0    0.4545
    1.0000         0    0.5455
    1.0000         0    0.6364
    1.0000         0    0.7273
    1.0000         0    0.8182
    1.0000         0    0.9091
    1.0000         0    1.0000];
colormap(map);

set(gcf, 'PaperOrientation', 'Portrait', 'PaperType', 'A4', 'PaperUnits', ...
    'centimeters', 'PaperPosition', [0.2 0.2 20.5 15]); 

load coast;

subplot('Position', [0.07 0.55 0.40 0.40]);
imagesc(phi, theta, CUr, [-.05 .05]);
hold on; plot([long; 360+long], 90-[lat; lat], '-k'); hold off; 
colorbar('h');
set(gca, 'Xtick', [0:60:360], 'YTick', [0:30:180]);
%title '(a) Re Hx; M2; sedi 5000';
title '(a) Re J_\phi';

subplot('Position', [0.07  0.05 0.40 0.40]);
imagesc(phi, theta, CUi, [-.05 .05]);
%imagesc(phi, theta, ((-Htr_1km_3 + Htr_1km_1)/max_Hxr_1km_1)*100, [-20 20]);
hold on; plot([long; 360+long], 90-[lat; lat], '-k'); hold off; 
colorbar('h');
set(gca, 'Xtick', [0:60:360], 'YTick', [0:30:180]);
title '(b) Im J_\phi';

subplot('Position', [0.56  0.55 0.40 0.40]);
imagesc(phi, theta, CVr, [-.05 .05]);
hold on; plot([long; 360+long], 90-[lat; lat], '-k'); hold off
colorbar('h');
set(gca, 'Xtick', [0:60:360], 'YTick', [0:30:180]);
title '(c) Re J_\theta';

subplot('Position', [0.56  0.05 0.40 0.40]);
imagesc(phi, theta, CVi, [-.05 .05]);
hold on; plot([long; 360+long], 90-[lat; lat], '-k'); hold off
colorbar('h');
set(gca, 'Xtick', [0:60:360], 'YTick', [0:30:180]);
title '(d) Im J_\theta';




print -dpng  'M2_Hx_1km_sedi_5000S_vs_no_sedi-sedi_5000S.png';

return